A theoretical analysis of flux-qubit measurements with a dc-SQUID 
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The readout process of a superconducting flux-qubit is theoretically analyzed in terms of the 
quantum dynamics of a qubit-SQUID coupled system during measurement. The quantity directly 
observed by the measurement is the switching current (7sw) of the dc-SQUID placed around the 
qubit ring. In order to clarify the relation between the /s„ and the qubit state, we calculated the 
time evolution of the density operator of the qubit-SQUID system while increasing the SQUID bias 
current until switching events occur and obtained the switching current distributions. This clarifies 
what information of the qubit is obtained by a dc-SQUID switching current measurement under 
specific conditions, for example, when the qubit eigenstate is a superposition of two different flux 
states. 

PACS numbers: 85.25.Dq, 03.67.Mn,03.65.Yz 



I. INTRODUCTION 



The quantum bit (qubit) is the fundamental ele- 
ment of quantum computergi. Recently, many experi- 
ments have been reported on quantum two-state systems 
consisting of superconducting circuits with Josephson 



junctions 



2.3.4.5.6.7.8.9.10 



Some of these systems use two 



flux-states in the superconducting ring that are macro- 
scopically distinguishable from each other, and the su- 
perposition of these two state3Sii^iSii£. These are called 
superconducting flux-qubits. The two states can be char- 
acterized by the directions of the circulating supercurrent 
along the ring, and the magnetic flux that is induced by 
the supercurrent. 

In experiments, small magnetic-field measurements 
made with dc-SQUIDs are used to read out the flux qubit 
states. A dc-SQUID is a highly sensitive magnetic flux 
probe. The measurement result is the bias current at 
which the SQUID switches to the voltage state. The 
switching current (Isw) of the SQUID varies depending 
on the quantum state of the qubit. 

The qubit ring has two stable states. One induces the 



flux 



and the other 



The SQUID switches at a 



bias current of /sw± = IcO cos[($sQ ± (t>q)/^o] depending 
on which state the ring is in. Here, <I>sq is the external 
flux applied to the SQUID ring, and $o = h/{2e) is the 
flux quantum. Therefore, at $sq/"I'o = V^: where the 
two states are energetically degenerated, the switching 
current appears probabilistically at /gw-i- or Isw- for each 
measurement in the sense of classical statistics. 

When Josephson junctions in the ring are small enough 
the ring becomes a quantum two-state system, that is, a 
qubit. Then, it is not so trivial what is meant by the 
obtained switching current, especially when the qubit is 
not in an eigenstate of the quantity measured by the 
SQUID. 

There are at least two possibilities. A flrst possibil- 
ity is that the dc-SQUID measurements of flux-qubits is 
a projection measurement of the smafl flux induced by 
the qubit ring current. The qubit has only two states 



where the current is definite. The switching current ap- 
pears probabilistically at /sw+ or Igw- even when the 
qubit state is a superposition of the two possible states 
before the measurement. Then, after the measurement, 
the qubit state jumps to the one of the two states corre- 
sponding to the measurement result. A second possibility 
is that the switching current has an intermediate value 
between /sw+ and Isw- and does not behave probabilis- 
tically even for a qubit that is in a superposition state. 



In experiments, when the fluctuation of the switching 
current is larger than |/sw+ — ^sw-| and when we need 
to take the average of many measured switching current 
values in order to determine the qubit-inducing flux, it 
is impossible to distinguish which of the two possibilities 
above is correct because the averaged values always ap- 
pear between /sw+ and /gw- independent of the possibil- 
ities and never behave probabilistically. To clarify what 
each switching current obtained by each measurement 
corresponds to, and what is meant by the distribution of 
a huge number of measured switching currents, we must 
undertake a quantitative analysis of the dynamics of the 
whole system, which consists of the measured object (flux 
qubit) and the measurement apparatus (dc-SQUID). 



The rest of this paper is organized as follows. Section 
II briefly introduces the background to superconducting 
flux-qubit measurement with a dc-SQUID. Section III 
shows the Hamiltonian of the flux-qubit and dc-SQUID 
coupled system, and derives a simplified version, which 
is applicable to numerical calculations. Section IV de- 
scribes our numerical calculation of the time evolution 
of the density operator of the qubit-SQUID composite 
system. Section V shows the results of the calculations 
and discusses them. Section VI makes some concluding 
remarks. 
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FIG. 1: Principle of the flux qubit. (a) Superconducting 
ring and fluxoid quantization, (b) Potential energy for the 
Josephson phase 7q and the two states, (c) Geometry of the 
dc-SQUID and the flux qubit for the measurement setup. 



II. FLUX-QUBIT AND READING OUT OF 
STATES WITH A DC-SQUID. 

A. Superconducting flux-qubits 

A superconducting flux qubit is a superconducting ring 
interrupted by a Josephson junction (or junctions). The 
fluxoid quantization causes the junction phase difference 
7q and the external flux $q applied by the magnetic field 
to obey the condition such that 



7q — 2mr + 2n 






(1) 



where $ = $q 



-^q-'cili 



n in an integer, and Idr is the 



circulating supercurrent along the ring [Fig. ^a)] . 

Suppose that an external flux $q ~ (n+l/2)<i>o pierces 
the qubit ring. The two states in the wells, <I>q ~ n$o and 
^q ~ {n+l)^o, of the potential energy F(7q) are energet- 
ically degenerate [Fig. ^b)]. We denote these two states 
by |L) and |R). The Josephson plasma oscillation energy 
hojp ~ y/2EcEj gives approximate level splitting in each 
well, where Ec = (2e)^/(2Cq) and Ej are the charging 
energy and the Josephson energy of the junction, respec- 
tively. When the potential barrier at ^^ ^ {n + l/2)<I>o 
is comparable to the charging energy Eq, macroscopic 
quantum tunneling of the phase couples states between 
the barrier. When the anti-crossing energy separation be- 
tween the coupled states are much smaller than the level 
splitting hujp, two states, those are energetically near- 
est and in different wells, are almost independent other 
states. Then, we can consider these two states as a two- 
state system, that is, a qubit. 

We can treat the states of the ring as a pseudo two-level 
system and obtain a reduced two-level-system Hamilto- 
nian with a {|L) , |R)} basis: 



i?o 



ea^ 



Act^ 



(2) 



where 2e is the energy difference between the localized 
states in the left and right wells and can be controlled by 



/. Here, A represents the transfer energy due to macro- 
scopic quantum tunneling between the wells, and 



a, = |L) {L\ - |R> {R\ 



a. = |L) {R.\ + |R) (L| ^ K 



1 
-1 

1 




(3) 



The eigenenergy of the ground (excited) states of 
Hamiltonian [3 is given by i?g(o) = 
eigenstates are 



-)V£2 + A2. The 



= sm 



|L> 



cos 



|L) -I- sin 



|R>, 



|R) 



(4) 



where tan[6'] = A/e. The energy levels exhibit anti- 
crossing at $q — {n + l/2)<I>o due to the mixing of 
|L) and |R) via macroscopic quantum tunneling. The 
quantum mechanical average of the qubit circulating cur- 
rent (/cir) = Ico{i\crz \i) [i = g, e), for the ground |g) 
and the first excited states |e), can be obtained using 
(g| (Jz |g) = — cos9 and (e| a^ |e) = cos0, where /co is the 
maximum supercurrent along the qubit ring. 



B. Measurement with dc-SQUID 

The "readout" of the flux qubit is the result of mea- 
surement to determine whether the qubit is in the |L) or 
|R) state. The difference between them, which is measur- 
able from the outside, is the direction of the circulating 
supercurrent along the ring and the magnetic field in- 
duced by the current. The field changes its sign depend- 
ing on the states |L) and |R). Then, we expect that, by 
measuring the field with a dc-SQUID, we can distinguish 
the state of the qubit. This method is actually used in 
the experiments. 

The switching current Isw of the SQUID is the bias 
current value at which the sum of the macroscopic phase 
differences of two junctions, 70+/ 2 = (71-1-72)72, escapes 
from the potential well. In a classical situation, where 
the external flux $sQ is added and a small flux ±0q 
is induced in the SQUID ring, the switching current is 
/sw(±(/)q) - /o(cos[7r$sQ/$o]T^sin[7r$sQ/*o]), where 
/o = (^) -Ejo, and i?jo is the Josephson couphng energy 
of the junction. In real qubit measurements, the induced 
flux is 10"3$o -^ 10"2$o. 

This type of quantum superposition state is called as 
macroscopic quantum coherence (MQC). The question of 
whether MQC is actually possible, was originally intro- 
duced by A. J. Leggett, in the 1980'si^. Recently MQC 
has been experimentally observed in a superconducting 
ring with Josephson junctions-. This MQC can also be 
used as a qubit because it is a typical two-state quantum 
system with controllable parameters. A qubit based on 
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FIG. 2: Quantum circuit for the flux-qubit and dc-SQUID 
coupled system. 



MQC in a superconducting ring is called a superconduct- 
ing flux-qubit. 

In experiments, the dc-SQUID is fabricated out of the 
qubit ring in order to measure a small magnetic field 
induced by the qubit circulating current [Fig. ^^c)] . The 
external field is tuned so that the fiux piercing the qubit 
was set near the degenerate point $q = (n + l/2)<I>o. 



III. DERIVATION OF SIMPLIFIED 
HAMILTONIAN FOR THE CALCULATIONS 

A. Total Hamiltonian for the qubit and dc-SQUID 
coupled system 

The real flux-qubit used in our experiments has three 
Josephson junctions. The dc-SQUID with two Josephson 
junctions has two leads. The ring and the dc-SQUID are 
coupled magnetically via mutual inductance M. This 
structure was proposed by Mooij 's group at the Technical 
University of DelfUi. 

The superconducting circuit for the total system that 
consists of the flux-qubit ring and the dc-SQUID is shown 
in Fig. 12 The Hamiltonian is given by 



H = H; 



SQ 



Ha + H, 



int 1 



(5) 



where 



Hn 



Co 



(7, 



-7-1- 2a7, 



q?) 



'2Ej cos[7q+/2] cos[7q_/2] - aEj cos[7q3] 

?iVL(7q-+7q3-2^/q)2 



+ 



2e 



2(LqL-M2) 



(6) 



is the qubit Hamiltonian, 



HsQ-^{^] ill^+lL) 



h 

2e, 

-2£;joCos[7o+/2]cos[7o-/2] - f^J ^bW^ 

^\ iq(7o- -27r/sQ)^ 



2e 



2{LqL- M^) 



(7) 



is the SQUID Hamiltonian, where Ih{t) is the bias cur- 
rent, which increases in time from zero to above the 
switching slowly enough compared with other relevant 
time scales, and 



Hu 



M (7o_ - 2^/sq) (7q- + 7q3 - 27r/q) 

(iqL - M2) 



(8) 

is the interaction between the qubit and the SQUID. 
Here, notations such as, 

7q± = 7ql ± 7q2, 7o± = 7oi ± 702 , (9) 

/q = $q/$0, /sQ = $Sq/*0, (10) 

are adopted. Moreover, L and Lq are the self-inductance 
of the SQUID ring and the qubit ring, respectively. M 
is the mutual inductance between them. Cq and Ej are 
the capacitance and the Josephson energy of two of three 
junctions in the qubit, and the other junction has values 
of aCq and ai?j with a ~ 0.8. Cq and Ejq are the ca- 
pacitance and the Josephson energy of the two junctions 
in the SQUID. 



B. Effective interaction between 70+, and cr 

1. Elimination of variables, 70- and 'y^s 

The Hamiltonian Eq. ((SJ of the total system contains 
five quantum variables; 70+, 70-, 7q-(-, 7q-, and 7q3. 
However, the quantity we really measure is the switching 
current /gw, that is, the SQUID bias current at which 70+ 
escapes from the Josephson potential well. Therefore, we 
have to analyze the way in which the relationship between 
the 7o+ of the SQUID and the qubit state evolves in time 
in the presence of the SQUID-qubit interaction. 

The qubit state is mainly determined by the value of 
7q_ , as described below. In terms of the qubit measure- 
ment, important variables are 7q_ and 70+. The other 
variables are not directly accessible quantities that can 
neither be controlled nor measured directly. Therefore, 
it is desirable that these variables are eliminated thus al- 
lowing the total system to be analyzed as a two-variable 
system. 

Fortunately, 70- and 7q3 can be approximated as har- 
monic oscillators strongly confined in parabolic poten- 
tials. We can integrate out these variables with the path- 
integral method, as described in Appendix A. Moreover, 



since the energy-level spacing of these harmonic oscilla- 
tors is very wide, the contributions, which are non-local 
in time, are negligible and we obtain a simplified Hamil- 
tonian for the total system as 



H' 



K 



^SQ 



H:„ 



(11) 



the ground state \g) and the first excited state |e) can be 
expressed as 



(18) 



|5)o - (|L) + |R»/V2, 
|e)o = (|L) - |R))/V2, 



where 



HL 



C„ 



7q++7q-+2a72_) 

2£'j cos[7q+/2] cos[7q_/2] - aE'j cos[7q_ - 27r/q], 

(12) 



where |L) and |R) are wavefunctions approximately local- 
ized in the left and right wells, respectively. When other 
excited states are adequately higher than these states in 
terms of energy, this is a two-state system. 

For finite but very small /, we can carry out a second- 
order perturbation calculation by adopting Hi as a per- 
turbation because f = fq — 1/2 is very small. Then we 
obtain the perturbed ground and excited states by the 
first order of perturbations as 



^SQ - ( ^ 



Co .2 

— 704 



- 2Ejo cos[7o+/2] cos[7r/sQ] 



Iff 



= l5)o 
^27TfaE, (^ 



(e| sin[7q_]|g)o 



e>o + E. 



o(»|sin[7q_]|g)o 
Ei 



i)o 



'2TrfaE3 



(- 



(g|sin[7q-]|e)o 
A 



l.9)o + E. 



o(»|sin[7q_]|g)o i 
Bi-2A I 



■f^int - 



1 

'2r„ 



(19) 



X |(iqa£;j sin[7q_ - 27r/q] 

-(Lqi - Al'')iEjo cos[7o+/2] sin[^/sQ])'} 



Here we disregarded Ejo, -Ej 

(?i/(2e))Viq, {n/{2e)f/L because 

{L and iq) are small in experimental situations 



where \i) and Ei is a second or higher excited state and 
MSjoCos[7o+/2]sin[7r/sQ])^i^^ energy when there is no perturbation. The qubit is 

always operated under the condition that the second and 
n4'\higher excited states are energetically far from \g)Q and 
|e)o, \Ei\^ \Ei — 2A| ;» A are always satisfied. Therefore, 
the third terms on the right hand side of Eq. (|19|) and 
the second and higher order contributions can be ignored. 
Using Eq. H18|l . we obtain 



compared with 
the inductances 



2. Two-state approximation of the qubit 

To achieve further simplification we have to calculate 
the quantum states of the three junction qubit using the 
Hamiltonian Eq. (|12|l with realistic parameters. In the 
vicinity of /q — 1/2, Eq. ((T^ becomes 



V2\^ 



|R) 



(20) 



where 



_ 2^(/q - l/2)aE, 

A 
2^(/q - l/2)aEj 



L|sin[7,-]|L) 



^q — -^0 



Hi, 



(15) 



where 



A 

x//^£(7g+,7g- 
^ £(/q)/2 



)V'L(7g+, 7g-) sin[7g_]d75+d75 



(21) 



Ho 



Ca 



(7q\ 



-^l_+2a^l_) 



and 



-2Ej cos[7q+/2] cos[7q_/2] + aEj cos[7q 



Hi = a£:jsin[7q_]27r/. 



(16) 



(17) 



and V'L(7g-i->7g-) = (7g-i->7g-|L)- Equation ^^ clearly 
shows that the ground and the first excited states can be 
expressed as superposition of |L)and |R) even when the 
external flux deviates slightly from the degenerate point 

/q = l/2. 

Therefore, we can simplify the qubit Hamiltonian H' 
as 



H' = £(/q)a, - Aa, 



(22) 



with / = /q — 1/2. Here, we ignored /^ and higher powers 
because /q is always in the vicinity of 1/2 (0.499 < /q < 
0.501) in flux-qubit experiments. Since Hq has an even 
symmetry about 7q_ and has potential wells on the both 
sides of 7q_ = 0, when Hi = 0, that is, when / = 0, 



where Uz and a^ are the spin operators with the basis 
{|L), |R)}. Although the above relations, Eqs. (|^ and 
()22|l . were obtained perturbatively, their validity is con- 
firmed by numerical calculations when e is comparable 

to A. 



This simplification is not solely for calculation con- 
venience; it also has a significant meaning in relation 
to the measurement process. Although in the original 
interaction Hamiltonian of Eq. (|14|) it appears that 
the SQUID measures the phase 7q_, on the real en- 
ergy scale, only the information about the qubit, |L) or 
|R), is transmitted to the SQUID. Therefore, it is now 
clear that the distribution of the switching currents ob- 
tained by enormous numbers of measurements never re- 
flects the probability distribution of the qubit wavefunc- 

tiOn|V'(7q-)P=l(7q-|V'q-)P. 

For further calculations, the interaction Hamiltonian 
Eq. H14|l should be expressed by this two-state approx- 
imation for the qubit. The qubit variable appearing in 
Eq. 1)14(1 is sin[7q_] and this should be expressed by the 
spin operators. Now, suppose that 



unitary evolution including the interaction in Eq. I|24|) . 
Since the bias current /b(i) changes with time, even 
the unitary evolution is calculated numerically with the 
forward Euler method. Although the 70+ value of the 
SQUID is a continuous variable, it is expressed by a dis- 
crete variable representation with a 31 ~ 41 basis. The 
prepared bases are the 70+ eigenstates when restricting 
the kinetic energy and the range of 70+ . 

Every pure state of the composite system can be ex- 
pressed as 



\i^^ 



d7o+ {0Lj(7o+)l7o+) ® |L) + 0Ri(7o+)|7o+) ® |R)} 

(26) 



where I70+) is the 70+ eigenstate, and (/)i(7o+) and 
4>r{io+) are the coefficients including the amplitude (i.e., 
not normalized). A mixed state is expressed as 



(R|sin[7q_]|R) = -(L|sin[7q_]|L)=7b, 
(L|sin[7q_]|R) = (R|sin[7q_]|L) = 0, 



(23) 



p{t) = X!p*i^*)^^» 



(27) 



are given. This shows that we can rewrite sin[7q_] — > 
7bfz- 7b is a constant of the order of unity, which is de- 
termined by the a and /q values of the qubit parameters. 
Finally, we obtain a simplified Hamiltonian for the 
qubit- SQUID coupled system 



where pi is the eigenvalues for the eigenstate \ipi) of the 
mixed state density operator. 

In order to show the switching current distribution, we 
plot the time evolution of the distribution function of 



H = e{f^)<Tz - Acr^j 



2EjQ cos[7o+/2] cos[7r/sQ] 



where 



m^-f 



1 
"2L0 



'2e\ ( 
— j |(Lqa£;j7b(T^ + MEjo cos[7o+/2] sin[7r/sQ]) 

-(LqL - Af2)(£;jo cos[7o+/2] sin[^/sQ])'}. 



p(7o+) = 

PL (704 
PR (704 



(7o+|Tr^[p(t)]|7o+), 
= Pl(7o+) +Pr(7o+), 



(28) 



Z^ 



Pi\(f>hiilQ+)\" 



^Pi\<l>Ri{H 



0+) 



IV. TIME EVOLUTION OF THE DENSITY 
OPERATOR p{t) 

A. Density operator for the qubit and dc-SQUID 
coupled system 

The simplified Hamiltonian for the qubit-SQUID cou- 
pled system contains two quantum variables, the qubit 
variable a and the SQUID variable 70+. Therefore, the 
density operator of the coupled system is defined with 
the basis of I70+, cr) and is expressed asi^ 

p{t) ^ d'yo+ dl'o+ Yl P'^'^'7o+,7i+Wl7o+,cr)(7o+, 

"' ■' (T,o-'=L.R 

(25) 
We use Feynman- Vernon's real-time path-integral 
methodic. We calculated the time evolution of the den- 
sity operator of the 7o+-o' composite system Using the 



(f>Ulo+) = ((7o+|«)(L|)|V'»), 

fe(7o+) = ((7o+|®(R|)IV'.)- (29) 

^ and Trg.^ [• • •] means taking the trace over the qubit 
states to obtain a reduced density operator for 70+. 
The switching current distribution Psw(-^b), which corre- 
sponds to the experimentally obtained distribution can 
be calculated from the portion of the SQUID wave- 
function remaining in the potential well pin(Jb(t)) = 
/i„d7o+(7o+|TrCT[p(i)]|7o+)by utilizing the following re- 
lation 



PsAh 



~A 



d/h 



^-A 



dpi, 



di 



W/d/b(i) 



di 



(30) 



where A is a normalization constant. Here, the region 
"in the well" is defined as — I.Stt < 70-1- /2 < tt. There- 
fore, J. d7o+ actually means the integration over the re- 
gion. The upper limit 70-1-/2 = n was chosen because 
.the energy barrier for the escape of the phase 70+ from 
the Josephson potential always has its peak value in the 
vicinity of 70-1-/2 = tt regardless of the amplitude of bias 
current. The speed d/b(t)/dt of the increase in the bias 
current is set at a constant value in examples discussed 
in the next section. 



B. Effect of decoherence 

Since we calculate the time evolution of the density 
operator of the qubit-SQUID coupled system, the effect 
of some types of decoherence can be easily taken into 
account. Here, we use a simple model of the decoherence, 
so called, the system-Boson modeliSii^. The Hamiltonian 
for the system and the Boson environment is given by 

Htot ^ H + HB+gAsJ2J2^bl + bk) + Y,^kblbk, (31) 

with the interaction Hamiltonian Hsb = Q^fzX]fc(^I + 
bk), where H is the Hamiltonian of the qubit-SQUID 
composite system given in Eq. (|24ll . bk is the annihila- 
tion operator of the Boson for the mode fc, and g is the 
coupling constant, and ^s is the operator of a physical 
observable in the qubit-SQUID coupled system, which in- 
teracts with the environment. The time evolution of the 
density operator expressed by the influence functional 
manner, and it can be calculated numerically by, for 
example, the program established by Makri's group-^. 
However, when it is allowed to approximate the effect 
of the interaction with the environment as a Markovian 
process, a simpler calculation can be performed, with 
the Lindblad equation^^. Although if the effect of the 
non-Markovian character of the system-environment in- 
teraction is important, the Lindblad calculation might 
not be very usefuUSiiS, however, it gives suggestive re- 
sults, as shown later, when the details of the origins of 
the decoherence in the experiments are unknown, as in 
the present case. 

The Lindblad equation is given by 

+f {2Asp{t)As - pit)AsAs - AsAsp(t)) , (32) 

where F is the coupling strength which is proportional to 
the square of g. 

When (i) As oc az, the decoherence destroys the quan- 
tum superposition between the \L) and \R) qubit states. 
This situation corresponds to the decoherence caused by 
the fluctuation of the external magnetic flux $q = /q^o- 
(ii) As oc 7o+ corresponds to the influence of the bias 
current /b the dissipation at normal resistance outside 
the SQUID. Moreover, (iii) As a d^o+/dt, expresses the 
detection of the finite voltage that appears across the 
SQUID. We will show the influence of this decoherence 
on the qubit-SQUID coupled system and the resulting 
switching currents by numerical calculations. 



V. RESULTS AND DISCUSSION 

A. Time evolution of the qubit-SQUID composite 
system during the measurement 

We carried out numerical calculations of the time evo- 
lution of the density operator given by Eq. (|25|l of the 



qubit-SQUID coupled system during the measurement. 

First, in Fig. |21 we show the time evolution without 
decoherence, in other words, the pure unitary evolution 
of the 7o+-cr2 composite system. Time t is normalized 
by the inverse of the Josephson energy Ejq. The pa- 
rameters for the numerical calculations are as follows. 
The Josephson energy and the capacitance of a SQUID 
junction are Ejq = lOOGHz, and Co = 40 fF, respec- 
tively. This gives Eq/Ejq ~ 0.01. The inductance of the 
SQUID ring L = 10 pH. The magnetic flux piercing the 
SQUID ring is $SQ = /sQ^o with /sq = 0.4. For the 
qubit, values of e = — 0.0065i?jo, and A — O.Oli^jo are 
given. These are almost the same as the experimental 
parameters used in the experiment by an NTT groupie. 
We employed the mutual inductance 7bAf = 1.4L, which 
gives approximately ten times greater qubit-SQUID cou- 
pling than in actual experiments carried out with the 
geometry shown in Fig. |21 Such strong coupling is cho- 
sen in order to distinguish the effect of the qubit state on 
the SQUID switching current within the precision of our 
numerical calculations. This interaction is still weak and 
does not change the discussion described later. 

This set of parameters corresponds to a situation where 
|L) and |R) are almost degenerate, (e 7^ because the 
qubit-SQUID coupling shifts the degeneracy point.) At 
t = 0, we set the SQUID bias current /b ^ 0, and the 70+- 
Uz composite system is in its ground state, that is, the 
bonding state of |L) and |R). For t > 0, /b is increased 
much more slowly than the time scale related to other 
relevant energies. Actually, in this and later examples, 
the bias current increases at approximately 400 nA/(/x 
sec). This corresponds to the bias current reaching the 
switching current from zero in a time of about 0.1 p, sec. 

The thin solid curves in Fig. |31 show the Josephson 
potential for 70+ at each time. The thick solid curves 
show PL (70+ ) , and the thick dashed curves show pr (70+ ) , 
which are defined in Eq. H29|) . The initial state has am- 
plitudes comparable for both the solid and dashed curves. 
This means that the initial qubit state is a superposition 
of |L)and|R). At this time, since 0l(7o+) oc 0r(7o+) is 
approximately satisfied, the qubit and the SQUID are not 
entangled. This is always the case when /b = regardless 
of the values of the parameters chosen. 

With increasing /b, the difference between |L) and |R) 
affects the position of the potential minimum in the co- 
ordinate 70-1-. As shown in the figure, without deco- 
herence the relation (^l(7o+) oc (/)r(7o+) is still largely 
maintained. Figure0]shows this fact quantitatively. The 
quantity 
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|/0£(7o+)0R(7o+)d7o4 



/d7o+|0L(7o+)| /d7o+|0R(7o+)r 



(33) 



for a pure state measures the degree of entanglement be- 
tween the qubit and the SQUID. R[^l,'4'r\ — 1, and 
indicate no entanglement (separable) and maximum en- 
tanglement, respectively. 

We can see that the quantity remains almost unity 
as the SQUID bias current /b increases. By using this 
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FIG. 3: Time evolution of the SQUID probability distri- 
bution without decoherence. A = 0.01-Ej. (a) Changes in 
the shape of the probability distribution. The solid curves 
show PL (70+) and the dashed ones show pr(7o+) defined by 
Eq. I29II . (b) Switching current distribution calculated by 
Eq. H3()|l . The time t is proportional to the bias current 
Ih{t), and normalized by the SQUID Josephson energy h/Ejo. 
d/b/di = 400 nA/(/i sec). 



quantity, the so called "entanglement of formation" is 
expressed as 

E = S{psq) = S{pquhit) = 

2 / d7o+|(/>L(7o+)|^ / d7o+|(?!)R(7o+)|^ (1 - i? [^l, 0r]) , 

(34) 

where pgq = Tr^[p] and Pqubit = Trsqip], and 
S{p') = Tr[p'logp'] is the von Neuman entropy. Since 
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FIG. 4: Time evolution of the measure R [<j!>l, ^r] of the en- 
tanglement between the qubit and the SQUID. R — 1 means 
no entanglement. 



2/d7o+|0L(7o+)l'/d7o+|'/'R(7o+)|' < 1/2 in Figs. 
and 21 the entanglement of formation E is always kept 
below 0.1%. 

The correlation between az and 70+ is formed by the 
interaction energy between the qubit and 70+. However, 
the entanglement makes it difficult to utilize the transfer 
energy A between |L) and |R) of the qubit because the 
SQUID wavefunctions (?!)l(7o+) and 0r(7o+) are different 
when entangled, and the overlap between them becomes 
small. It is postulated that this fact suppresses the en- 
tanglement formation. In this situation the switching 
current may have an intermediate value between the val- 
ues for |L) and |R). The switching current distribution 
is shown in Fig. |2l[b). It has only one significant peak. 



B. What do we measure Avith a dc-SQUID ? 

In the vicinity of a half integer of /q, the ground state 
\g) and the first excited state |e) of the qubit are 



Iff) = a(/q)|L) + 6(/q)|R), |e) ^ 5(/q)|L) - a(/q)|R), 



respectively. Here, we define 



(35) 



a(/q) 



(/q) 



D 



HU) 



D 



(^(/q) + \/ 



e(/q)2 + A2 



(36) 



Since e(/q) = at /q = 1/2, \g) and |e) are superposi- 
tions of |L) and |R) at /q ~ 1/2. Therefore, a^, which is 
proportional to the flux induced by the qubit circulating 
current, has large quantum fluctuations, and does not 
have a definite value. Then, what the dc-SQUID switch- 
ing current Isw measuring the induced fiux corresponds 
to is ambiguous. 

A simple explanation of the ct^ projection measure- 
ment is as follows. It requires not only classical correla- 
tion but also quantum entanglement between the SQUID 
and the qubit. When the interaction between the qubit 
and the SQUID produces entanglement between them, 
as given by Eq. (|26|l . the wavefunctions (/)Li(7o-i-) and 
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FIG. 5: Schematic of dc-SQUID switching current distribu- 
tion according to a "projection" measurement of the qubit. 
(a)Before averaging. (b)Average value. 



0Ri(7o+) have their peaks at different positions, 70L and 
7oR. As the SQUID bias current is increased, and when 
one of their wavefunctions reaches its switching point, it 
switches probabihstically. Therefore, if the two peaks are 
adequately separated in the coordinate 70, we may find 
two peaks in the switching current distribution. 

We denote the /sw for |L) and |R) as /l and /r, respec- 
tively. According to a simple measurement formulation, 
the above is a a^ measurement. Then, if the simple for- 
mulation were correct, the measurement should give 

• /l with the probability |a(/q)p, /r with the prob- 
ability |6(/q)|2. 

This is a probabilistic phenomenon, and therefore, the 
obtained values have a large distribution. These two 
peaks positioned at /l and /r appear in the plot of a huge 
number of obtained values. After averaging the values, 
the expectation value 

(/sw) = (5|/sw|g) = |a(/q)|'/L + |&(/q)l'/R (37) 

would be obtained. This is shown schematically in Fig. 

El 

If the measurement with the SQUID were a projec- 
tion measurement of g^ as described above, the observed 
switching current would appear at a value correspond- 
ing to 0l(7o+) or (/)r(7o+), probabilistically. There- 
fore, if 0l(7o+) were not proportional to 0r(7o+) and 
i?[0L,0R] ^ Ij the switching current would split into 
two values even for a one qubit state at absolute zero 
temperature, as sketched in Fig. [SJa). 

In the experiment reported by the Delft group^, (per- 
haps for a reason that has nothing to do with qubit quan- 
tum fluctuations) the /gw had a broad distribution, which 
prevents one even from distinguishing /l and /r on a sin- 
gle measurement basis. Therefore, they showed only the 
averaged values, which faithfully reproduced the curve 
shown in Fig. EJb). 

However, as described above, when we carried out nu- 
merical calculations on the time evolution of the density 
operator of the total system, we found that the total sys- 
tem wavefunction with low energies does not adopt an 
entangled state but remains almost separable, that is, 
-R [0L, -^r] ~ 1 (</'l(7o+) oc 0fl(7o+) = </'(7o+))- Only the 
position in the 70+ coordinate of the wavepacket 0(70+) 
depends on the qubit state. The obtained switching 



current corresponds to that of (?!)(7o+), which is single- 
valued, unlike the entangled state mentioned above. In 
addition, the numerical calculation also showed that the 
switching current has a two- value distribution in the pres- 
ence of strong decoherence, which destroys the qubit su- 
perposition, as we will discussed in Figs. |3 and ^2 
Since the qubit circulating current /dr is proportional 
to the z-component of the qubit spin, it should be writ- 
ten as /cir = /coCzj under the two-state approximation 
for the qubit. According to quantum mechanics, when 
the spin is in a superposition state |'0) = a |L) + 6|R), 
each Oz measurement provides the discrete result -1 or 
1 probabilistically, and never provides an intermediate 
value. Although the quantum-mechanically averaged 
value of the circulating current for the superposition state 

1^) = a |L)-h5 |R) is given by (Jdr) = /co (2 \af - l) , this 

value is obtained after we measure lea for the state |V') 
many times and calculate its average value. 

Very recently, the NTT group has succeeded in per- 
forming single-shot detection of the qubit flux by improv- 
ing the switching current resolution in their SQUID —. 
The plots of the raw values of the switching currents 
clearly show the qubit state transition from the ground 
state to the first excited state accompanied by a change in 
the external magnetic field <I>q = fq^o, without any av- 
eraging of the dataiS. One of the most interesting facts 
in their data is that the switching currents illustrate a 
curve that explicitly shows that the qubit is in a super- 
position of |L) and |R) in the vicinity of a degenerate 
point /q — 1.5. In such a region, their switching cur- 
rents appear at intermediate values between the switch- 
ing currents for |L) and |R). A curve that is very similar 
to that in Fig. ]Ej[h) appeared from the single shot data. 
The behavior of the switching current versus the exter- 
nal flux / in Fig. |3Jb) is very similar to that of the 
quantum-mechanically averaged value of the qubit circu- 
lating current (/cir) in Fig. 2(b). At first glance, this is 
very strange if we adopt the above simple formulation of 
the dc-SQUID measurement, which gives a two peak Igw 
distribution for a superposed qubit state. 

However, taking an experimentally realistic weak 
qubit-SQUID interaction into account in the Hamilto- 
nian of the composite system, one finds that no entan- 
gled state between the qubit and the SQUID is formed. 
The time evolution keeps the composite system to be in 
a direct product state even in the presence of the inter- 
action; 



IV'tot) = IV'Sq) «) \Tpq 



(38) 



The qubit state affects the shape of the SQUID wave- 
function; 



H10+) = (to+I-'/'sq)- 



(39) 



The position of the peak of the wavefunction 0(70+) 
varies with the qubit state although the wavefunction 
is a wavepacket with a finite width. 



C. Variational analysis of the entanglement 
formation between the qubit and the SQUID 

Entanglement formation in time is too complicated to 
examine analytically. Therefore, here, we consider the 
behavior of the entanglement in the ground state of the 
qubit-SQUID coupled system by using a variational cal- 
culation. 

We adopt the trial wavcfunction 

IV'tri) - — / ( 6Le ^ |7o+) ® |L) 



+6Re 



2 
7ra/ 

(70+/ 



-mr 



1 70+; 



|R) d7, 



o+, 



(40) 



for the total system because the SQUID potential is 
nearly parabolic and the wavcfunction in the potential 
should take the ground level. The width of the wave 
packet is a. We denote the difference between the peak 
positions of the L and R components as d = {7L — 7r)/2, 
and mid point between them by 7 = (7L + 7r)/2. From 
the Hamiltonian (|24|l . the expectation value (V'tril^lV'tri) 
of the energy for the trial wavcfunction is given by 

V 

(V'tril-fflV'tri) = 

- 2e^°'I^EiQ cos[7r/sQ] {cos 70 cos d 
-(|6l|^- |feR|^)sin7osin(i} 



_c 
a 



-a/8, 



£;jofc{sin7osind- (I^lI^ - |6rP) cos7n cosd} 
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+ e{H\'-\bn\'), 



where 



Ma f2e\ , ^ , 

k = -^ I y ) ^J7b sm[7r/sQj, 



(41) 



(42) 



is the dimensionless qubit-SQUID coupling constant, and 

Here k is of the order of 10~^ when we use the parameters 
in experimentsiiS. Therefore, we disregard the higher 
powers of k in the calculations below. Here, we ignored 
the self-energies for the qubit and the SQUID induced by 
the interaction, to simplify the discussion of the entan- 
glement formation. 

We minimize Eq. (|41|l with respect to the peak posi- 
tions 7l, and 7r,, and the width of the wavcfunction a, 
using an approximation that disregards k^, cP, a?' and 
higher powers. Then, we obtain the minimum 
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2 cos[7r/sQ] 



(45) 



(46) 



(47) 



This minimum is realized when the peaks are positioned 
at 



7o+ 
2 



= 7l = 7o 



da, — =7r = 7o -ao, 



(48) 



where 
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(51) 



This split between peaks do (A 7^ 0) is much smaller than 
the difference between the peak positions do (A — 0) for 
the simple (not superposed) |L) and |R). In fact, the 
ratio is given by 



do 



do(A = 0) 



VT^^Vh^ - h^ 



^/T^^yV^V + ^(f) 



(52) 



This ratio becomes vanishingly small as the bias cur- 
rent approaching the switching current, that is, /b — > /o- 
Therefore, when the qubit superposition A exists, the 
entanglement between the qubit and the SQUID is neg- 
ligible, at least just before the switching. At that time, 
Eq. (|49|l gives the position of the SQUID wavcfunction 
peak at 



70 



arcsm 



/h 



(1 - k'x) 



(53) 



^"j ({I + xk')^Io'' - hit)^ + /b(i)arcsin[/b(i)//o]') 



This coincides with the peak position when the qubit has 
a classically definite spin az = x. 

Finally, we can estimate by an approximate variational 
calculation that the peak of the SQUID wavcfunction is 
positioned at 
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FIG. 6: Schematic dc-SQUID switching current distribution 
according to a "weak" measurement of the qubit. (a)Before 
averaging. (b)Average value. 



when the qubit is in a superposed state |V'q) = a|L) +6|R) 
unless the state of the qubit undergoes a rapid change. 
This coincides with the quantum mechanical average 
(V'ql^zlV'q) of 7o+- Therefore, in the absence of deco- 
herence, we can interpret this to mean that the resulting 
switching current value, at which the 70+ escapes from 
the Josephson potential, coincides with the quantum- 
mechanical average (/sw) before any averaging operation. 
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D. Effect of decoherence 



L=91520, -1^0.239556 



The discussion in Sec. IV Bl becomes more convincing 
when we compare the time evolution in Fig. 0with de- 
coherence. 

When we use unrealistically strong coupling between 
the qubit and the SQUID, there is a significant entan- 
glement. However, this still gives only one peak in the 
switching current distribution in our calculation. This 
shows that entanglement by itself does not cause a two- 
value measurement, in other words, a a^ projection mea- 
surement. Nevertheless, such strongly entangled states 
are not robust against a weak decoherence, which de- 
stroys the qubit superposition and realistic switching cur- 
rent distribution for a strong interaction as discussed in 
this section. 

Assuming a Markovian process, we performed the cal- 
culation using the Lindblad equation (|32|l with As = (Jz 
and r = 3.6 x 10~^£'jo/?i. The decoherence makes the 
composite system a mixed state, and the states |L) and 
|R) become incoherent to each other. Therefore, the 
transfer energy between |L) and |R) through A does not 
matter any longer, and the interaction between the qubit 
and the SQUID determines the state of the system result- 
ing in the qubit state dependence of the SQUID wave- 
function shape. In this situation, however, the qubit- 
SQUID correlation is not an entanglement but a classical 
correlation because the states with different qubit states 
are already incoherent to each other. In terms of Eqs. 
(|77|l - (|^ of the density operator, the L component cpi^i 
and the R component (/)pjj belong to different \ipi), IV'i) 
{i ^ j). The switching current value becomes one for |L) 
or for |R) in a probabilistic manner. This appears in the 
switching current distribution in Fig. [7fb). The peak is 
split into two, corresponding to |L) and |R), respectively. 

This type of peak splitting in the switching current 
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FIG. 7: Time evolution of the SQUID probability distri- 
bution with decoherence proportional to a^ with F — 3.6 x 
lO^^'^Ejo/h. A = 0.001-Ej. (a) Changes in the shape of prob- 
ability distribution. The solid curves show pl(7o+) and the 
dashed ones show pr(7o+) defined by Eq. 12911 . (b) Switching 
current distribution calculated using Eq. H3()|l . The time t is 
proportional to the bias current /b, and is normalized by the 
SQUID Josephson energy Ti/ Em. d/b/df = 400 nA/{fi sec). 



distribution caused by the decoherence occurs when the 
qubit superposition is weak (A is small), or/and when 
the qubit-SQUID interaction energy is so large that a 
strong qubit-SQUID entanglement is formed before the 
decoherence makes the system a mixed state. 

Our formulation can readily be used to investigate an- 
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other type of decoherence. Although we do not show the 
results of our numerical calculations here, we briefly com- 
ment on some of them. When the quantity of the system 
that couples with the environment is the SQUID 70+ or 
7o+ , the interaction between the system and the environ- 
ment does not directly deteriorate the qubit superposi- 
tion between |L) and |R). Since the qubit-SQUID cou- 
pling through the mutual inductance is weak in the above 
calculations (and experiments), the environment deforms 
only the shapes of the SQUID wavefunction. This makes 
the qubit-SQUID system a mixed state. However, each 
eigenstate of the mixed-state density operator of the 
|L) — |R) superposed qubit is largely maintained even if 
the coupling to the bath is strong. Although the switch- 
ing current distribution of the mixed-state density oper- 
ator may sometimes have several peaks, the splitting of 
the peak does not correspond to the splitting of the qubit 
into |L) and |R). The peaks instead correspond to dif- 
ferent SQUID states, for example, the ground and first 
excited states of the Josephson plasma oscillation. 



E. Switching current behavior with a change in the 
external magnetic flux cf>q 



If the above discussion of what is indicated about the 
qubit by the switching current is correct, our theory must 
reproduce the switching current behavior observed in ex- 
periments with the change of the external magnetic flux 
^q (oc /q) shown in Fig. |2| We carried out a time- 
evolution calculation at various /q values and present the 
results in Figs. 013 and CHI 

Figures IHl and shows the switching current behavior 
with the energy splitting e between |L) and |R) when 
there is no decoherence. Since e is proportional to the 
deviation of the external magnetic flux from the degen- 
erate point fq — {n + 1/2), as shown in Eq. H21|l . the 
horizontal axis also corresponds to the change in fq. The 
peak position of the switching current distribution grad- 
ually shifts from /l to /r, and there is no discontinuous 
transfer. This suggests that the switching current traces 
a continuous curve as shown by the schematic illustration 
in Fig. m 

Figure ^2 where the peak positions of the switching 
current distributions for various external fluxes are plot- 
ted, shows clearly that the e dependence of Igv, obeys 
such a law. Without decoherence Igw gradually transfers 
from the switching current for |L) to |R) across the de- 
generate point (e = 0). By contrast, Igw suddenly jumps 
from one for |L) to |R) at the degenerate point and two 
peaks appear at certain flux values. The circles indicated 
by arrows are the Igw values of the same parameter sys- 
tem for various decoherence strengths. This provides ev- 
idence that a measurement giving a single intermediate 
Isw is changed into a two-value measurement by strong 
decoherence. 
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FIG. 8: Calculated dc-SQUID switching current distri- 
bution in the qubit measurement. (A — O.Ol-Bjo) Each 
curve corresponds to the distribution for an energy dif- 
ference £ between |L) and |R). The curves are verti- 
cally shifted for clarity. From top to bottom, e/Ejo = 
-0.025, -0.015, -0.010, -0.0065, -0.0046, -0.002, -^0.005. 
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FIG. 9: Calculated dc-SQUID switching current distri- 
bution in the qubit measurement. (A — 0.01-Bjo) Each 
curve corresponds to the distribution for an energy dif- 
ference e between |L) and |R). The curves are verti- 
cally shifted for clarity. From top to bottom, e/Ejo = 
-0.027, -0.025, -0.015, -0.010, -0.0065, -0.004, -0.002, 
+ 0.005, -F0.007. 
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FIG. 10: Calculated dc-SQUID switching current distri- 
bution in the qubit measurement. (A = O.OOl-Ejo) Each 
curve corresponds to the distribution for an energy dif- 
ference e between |L) and |R). The curves are verti- 
cally shifted for clarity. From top to bottom, e/Ejo = 
-0.015, -0.010, -0.0065, -0.004, -0.002, 0.000, +0.002, 
+ 0.005. 
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FIG. 11: e dependence of the switching current peak position. 



F. Comment on the precision of a single 
measurement 

The problem analyzed here is invoked by our experi- 
ment where each single switching current measurement 
has sufficient resolution to distinguish whether the qubit 
is in |L) or |R). If we could not obtain the switching cur- 
rent behavior depending on the external flux, it would 
not make us aware of the importance of understanding 
the switching currents that appear in the intermediate 
region between /L(for |L)) and /R(for |R)) because the 
averaged value of the projection measurement results al- 
ways gives the intermediate value when the qubit is in a 
superposition state. Therefore, here, we briefly consider 
the factors that determine the resolution of the switching 
currents. 

As we have discussed in this paper, the quantum fluc- 
tuations in the qubit-induced flux (/)q, in other words, 
the quantum fluctuations of cFz , are not directly related 
to the SQUID resolution. Within the precision of our 
numerical calculation, there is no significant difference 
between the half width of the switching current distri- 
bution for a Uz definite state and a superposition state. 
In our numerical calculation, we found that the distri- 
bution width increased when we used a larger SQUID 
junction capacitance. When the capacitance is large, the 
level splitting of the SQUID plasma oscillation of 70+ is 
narrow, so the increase in the SQUID bias current ex- 
cites the plasma oscillations. As the result, the state 
of 7o+ is distributed over many levels that have slightly 
different switching currents. Therefore, the distribution 
becomes broader resulting in a poorer resolution. This 
fact does not identify the origin of the high resolution of 
our measurement because there are many other factors 
that cause those changes in the distribution width, in- 
cluding the speed of the increase in the bias current, and 
the strength of dccohcrcncc. The resolution problem is 
still an open question. 



VI. CONCLUSION 

We discussed the superconducting flux-qubit measure- 
ment with a dc-SQUID from the viewpoint of the dy- 
namics of a composite system consisting of a qubit and 
a SQUID. Since both the qubit and SQUID are quan- 



tum systems the relation between the measured object 
(qubit state) and the outcome (switching current) of the 
measurement apparatus is not trivial. 

We calculated the time evolution of the density oper- 
ator of a qubit-SQUID coupled system from the onset of 
an increase in the bias current to the switching event, 
and examined the behavior of the switching current dis- 
tribution. 

Of particular interest is the fact that when the ground 
state of the composite system is a qubit superposition 
state, the measurement designed to distinguish in which 
state the qubit is does not always work as designed, and 
there is a possibility that the switching current will have 
an intermediate value between the two states. Never- 
theless, even in such a situation, when there is decoher- 
ence (reduction in the off-diagonal elements of the density 
operator) which destroys the qubit superposition of |L) 
and |R), the switching current then gives us information 
about |L) or |R) in a probabilistic manner. 

When the external field $q is far from the degenerate 
point ($q = {n -\- l/2)$o), since there is no quantum 
transfer between the two states of the qubit, the qubit 
should be in the lower energy state, and the switching 
current of the SQUID shows a value corresponding to 
that state. The switching current Igw is different for 
<I>q < (n + l/2)<I>o and $q > {n + l/2)$o because the 
energies of the two states cross at the degenerate point. 

We find by variational calculation that the peak posi- 
tion of the wavepacket (/)(7o+) is approximately the bot- 
tom of the SQUID Josephson potential Vsq(7o4.;(Tz = 
2|ap — 1). The resulting switching current becomes that 
for a usual magnetic field measurement with SQUID of 
this potential. This value coincides with the quantum- 
mechanical average (/sw) of the switching current for the 
qubit j-j/iq) — a |L) -I- 6|R). The fact that the switching 
current behaves as the average value, which may at first 
seem strange, is explained as described above. Moreover, 
this indicates a very important fact about the dc-SQUID 
measurement of flux-qubit states. The measurement is 
not a projection measurement that determines whether 
the qubit is |L) or |R). It is a simple measurement of the 
weight |ap(= 1 — |6p) of the qubit. 

The flux 0q induced by the qubit ring current has a 
large quantum fluctuation when the qubit is in a super- 
position state. The fluctuation of the measured switching 
current, however, has almost nothing to do with the fluc- 
tuation of (/)q . The switching current fluctuation is caused 
mainly by fluctuations of the SQUID itself. 

Our theoretical results in the present work clearly ex- 
plain the switching current behavior revealed by our ex- 
perimental results. This may provide us with important 
clues to improve the measurement scheme for quantum 
computations. 
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APPENDIX A: ELIMINATION OF VARIABLES, 

70- AND 7q3 

Here, we briefly explain how we eliminate invisible vari- 
ables, 7o- and 7q3. 

Picking up terms including 70- or 70- from the Hamil- 



tonian Eq. Q, we obtain 

i?(7o_) = -TOi;^ + -mci;o^a;^-|-xF(i)-2^joCos[7o+]sin[7r/sQ] 



where 



Co fh 



2 V2e 



7o- - 27r/sQ, 



(Al) 
(A2) 
(A3) 



2Ln 



LUo 



Co(LqL - M2) 



2e\ 2E 



JO 



Co 



cos[7o+/2]cos[7r/sQ], 
(A4) 



(A5) 

7q = 7q- + 7q3 ~ 27r/q. (A6) 

Here, we neglected x'^ and higher powers such as cos[a;] — > 
1 — x^ jl, sin[a;] — > x because x is strongly constricted in 
the parabolic potential. 

Then, the Feynman- Vernon influence functional ex- 
pressing the effect of 70- on the motion of F(t) is given 
b>i^i^ 



where 



T^M = exp [-$Fv] , 



$FV = - / dti / dt2[F(ii)-F'(tl)] 



(A7) 



X [L(ti - i2)F(t2) - i*(ti - t2)F\t2)\, (A8) 



L(t)^ 



IniLOo 



coth 



hujof3 



cos[clIo^] ~ isiii[u!ot 



(A9) 

Carrying out the ^2 integration in Eq. IjASp by part 
two times, 



dt2 [Lit, - t2)Fit2) - L*it, - t2}F'{t2)] 



-iiF{ti) + F'(ii) - cosujoti{F{0) + F'(0)) 



2mujQ 
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+ 



coth 
1 



hujof] 



smu;oti{F{0)-F'{0)) 



2mwo 



-i sinLJoti{F{0) + F'(0)) + coth 



hujoP 



{coscooh{F{0) - F'{0)) - {F{h) - F'ih)) 



^0 



Lih - t2)F{t2)At2 



L*{h-t2)F'{t2)dh 



(AlO) 



is given. Here, F(t) indicates the quantity on the for- 
ward hne and F'{t) that on the backward one. 'A' 
and 'A' mean the first and the second time deriva- 
tives of the quantity A, respectively. Since F{t) = 

iTef-L^r^Ut) - ^sin[7o+/2]si n[7r/sQ ]7o+, and 
both 7q(t) and 70+ are of the order of ^/eJEc/H at most, 
the terms from the second line of Eq. IJA10|I are smaller 
by the factor \/EJE^ / (hujo) ~ (^) \/LEjo than those 
in the first line. This factor is no larger than 10~^ in 
experiments^i^iiS, therefore, we only take the terms in 
the first line in Eq. (|A10|) . 
Now we obtain 



$ 



FV 



2hmL0Q jQ 



dh[Fih)-F'{tl)] 



-iiFih) + F'ih) - cosuJohiFiO) + F'{0)) 



-coth 



suYUJoti{F{Q)-F'{Q)) 



(All) 



Highly oscillatory terms such as, sinwo^i, costJoii are 
suppressed by the integration over ii. This enables us to 
neglect these terms. 

Finally, we obtain the influence functional 



^¥Y 



w exp 
= exp 
= exp 
X exp 



2hmujQ Jo 



M7q(^) 
Lqi - Af 2 



Sjo COS [70+ /2] sin[7r/sQ] 



(A13) 



Very similarly, the elimination of 7q3 is carried out as 
follows. The Hamiltonian concerning 7q3 is 

■f^(7q3) = i^myv^ + -TOyWi^y^ + Fy{t)y 

2in 



^) ^jo'cos2[7o+/2]sin2[^/sQ] 



aE:^ cos[7q_ - 27r/q] - (7q_ - 27r/q) x 
a£:j sin[7q_ - 27r/ql -^ 



1 



ir 



aEi sin[7q_ - 27r/q] 



E'jo cos[7o+/2] sin[7r/sQ] 



1 
2L„ 



where 



^Cn [^ 



y = 7q3, 



lyUi = aEj sin[7q_ - 27r/q] + -;- { tt 



(7q--2^/q)2, 

(A14) 

(A15) 
(A16) 
, (A17) 



Fy{t) 



(7q- - 27r/q) I aE] cos[7q 
M 



2-/q]-(^) - 



dh[F{h) - F'{ti)][F{h) + F'ih)] 



Ejo cos[7o+(t)/2] sin[7r/sQ] - aEj sin[7q_ - 27r/q], 

(A18) 



2h'mujQ Jo 



2h'mujl Jo 



dhiFihf-F'ihY] 



F{tiYdU 



2fimuj^ Jo 



F\hYdh 



This influence functional is local in time and F{ti) and 
F'{ti) are separated so that the effect of 70- gives a 
unitary evolution which is expressed by the Hamiltonian 



Here, we neglected (7q3-|-7q_ — 27r/q)^ and higher powers 
because they are strongly constricted in the parabolic po- 
tential. The influence functional expressing the effect of 
7q3 on 7q_ and 704- is given from this description. More- 
over, the terms that are non-local in time are negligible 
, . for the same reason as that for 70- . Finally, we obtain 
tne simplified Hamiltonian Eq. (|12|l where 7q3 and 70- 
have been eliminated. 



H 



2mujQ 



F{tf 



L^L ~ AP 

2r„ 



APPENDIX B: DISCRETE VARIABLE 
REPRESENTATION OF 70+ 

Since the variable 70+ is a continuous one, it is not 
directly suitable for numerical calculations. Therefore, 
we discretized them as follows. 
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The domain of 70+ is restricted within (a, 6), where 
a = — 0.77r,6 — 2.57r. This domain is wide enough to ex- 
press the SQUID wavefunction. We prepare orthonormal 
basis {|fc„)}, which satisfies (7o+|fc„) = exp[ifc„(7o+ — a)], 
where fc„ = 2m: /{h - a), n = 0, ±1, ±2, • • • , ±iV. They 
are Fourier series and 70+ cigcnstatcs. From this set of 
{|fc„)}, we make 70+ cigenstates 

7o+l7m>=7m|7m>, (m - 0, ±1, ±2, • • • , ±iV), (Bl) 

where 

|7m) = ^flnmlfen), ^ \amn\^ = 1- (B2) 



We use this basis {|7m)} for the discrete variable rep- 
resentation of 7o+ and the wavefunction of the SQUID is 
expressed as the hnear combination of \^m)- Using this 
basis with finite N corresponds to the abandonment of 
the states of higher momentum energies in the SQUID. 
Therefore, we can calculate low energy states precisely 
with small numbers of elements of the basis compared 
with other bases, such as 70+ eigenstates placed with the 
same interval over (a, 6). 



